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ABSTRACT 

Gravitationally lensed quasars can be used as powerful cosmological and astrophysical probes. We can (i) infer the Hubble constant H 
based on the so-called time-delay technique, (ii) unveil substructures along the line-of-sight toward distant galaxies, and (iii) compare 
the shape and the slope of baryons and dark matter distributions in the inner regions of galaxies. To reach these goals, we need high- 
accuracy astrometry of the quasar images relative to the lensing galaxy and morphology measurements of the lens. In this work, we 
first present new astrometry for 11 lenses with measured time delays, namely, JVAS B02 18+357, SBS 0909+532, RX J091 1.4+0551, 
FBQS J0951+2635, HE 1104-1805, PG 1115+080, JVAS B1422+231, SBS 1520+530, CLASS B1600+434, CLASS B1608+656, 
and HE 2149-2745. These measurements proceed from the use of the Magain-Courbin-Sohy (MCS) deconvolution algorithm applied 
in an iterative way (ISMCS) to near-IR HST images. We obtain a typical astrometric accuracy of about 1-2.5 mas and an accurate 
shape measurement of the lens galaxy. Second, we combined these measurements with those of 14 other lensing systems, mostly from 
the COSMOGRAIL set of targets, to present new mass models of these lenses. The modeling of these 25 gravitational lenses led to the 
following results: 1) In four double-image quasars (HE0047-1746, J1226-006, SBS 1520+530, and HE 2149-2745), we show that the 
influence of the lens environment on the time delay can easily be quantified and modeled, hence putting these lenses with high priority 
for time-delay determination. 2) For quadruple-image quasars, the difficulty often encountered in reproducing the image positions to 
milli-arcsec accuracy (astrometric anomaly problem) is overcome by explicitly including the nearest visible galaxy/satellite in the 
lens model. However, one anomalous system (RXS Jl 131-1231) does not show any luminous perturber in its vicinity, and three others 
(WFI 2026-4536, WFI 2033-4723, and B2045+265) have problematic modeling. These four systems are the best candidates for a 
pertubation by a dark matter substructure along the line-of-sight. 3) We revisit the correlation between the position angle (PA) and 
ellipticity of the light and of the mass distribution in lensing galaxies. As in previous studies, we find a significant correlation between 
the PA of the light and of the mass distributions. However, in contrast with these same studies, we find that the ellipticity of the light 
and of the mass also correlate well, suggesting that the overall spatial distribution of matter is not very different from the baryon 
distribution in the inner ~5 kpc of lensing galaxies. This offers a new test for high-resolution hydrodynamical simulations. 

Key words. Gravitational lensing: strong - Galaxies: quasars: general - Techniques: image processing 
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^ . 1. Introduction J2009I I201O) showed that a detailed study of the lensed quasar 
5— i CLASS B1608+656, combining different high-resolution data 
d ; The measurement of the time delay between the lensed im- and advanced lens mo deling techniques, could lead to an esti- 
ages of a gravitationally lensed quasar offers one of themost mate of Hq with smal] random and systematic uncertain- 
elegant ways to measure the Hubble constant, H (|RefsdaJ| ^ Hq = 7Q 6 ± 3 j km s -i Mpc -i. The determination of H 
[1964D. This cosmological application of strongly lensed quasars based on a large sample of lensed quasars can also reach a pred . 
motivated many of the early gravitational lensing studies and sion competitive with the one of standard methods (0guri 20 07; 
time-delay measurement campaigns. Unfortunately, the preci- q q \ ss 2008) 

sion on Hq obtained with the time delay method has been, un- T , . . , , , , 

1 11 j 1 • r. 1 \ 1 1 In order to use the time-delay method as a precision- 

til now, challenged by systematic errors. Recently, Suvu et al. , , • • . • • ■ , 

J J J — cosmology tool, it is necessary to derive accurate time-delay 

* Based on observations made with the NASA/ESA HST Hubble measurements and to properly understand the uncertainties asso- 

Space Telescope by the CfA-Arizona Space Telescope Lens Survey clated t0 the lens environment and to the lens mass distribution. 

(CASTLeS) collaboration, obtained from the data archive at the Space The COSMOGRAIL project, i.e. the COSmological MOnitoring 

Science Institute, which is operated by AURA, the Association of of GRAvltational Lenses, aims to provide the community with 

Universities for Research in Astronomy, Inc., under NASA contract exquisite lightcurves and accurate time-d elay measurements for 

NAS-5-26555. a sample of more than 20 lensed quasars (Eigenbrod et alJ[2005b 
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IVuissoz etal J 120071 l2008t ICourbin et all 12010). To understand 
degeneracies affecting lens models, we started to model all sys- 
tems monitored by COSMOGRAIL as well as those monitored 
by earlier campaigns in a uniform way. Because the relative 
astrometry of the lensed images is the main constraint on the 
lens models, we first worked on deriving an updated astrometry 
of these systems by applying the Itera tive Strategy coupled with 
the MC jj deconvolution algorith m (Magain et al. 1998), i.e. 
ISMCS dChantrv & Magainl l2007) to HST images. This allows 
us to obtain accurate relative astrometry and morphological 
information of the lensing galaxy (i.e. ellipticity and position 
angle PA of its major axi s ). We presented the first part of 
this work in Chantrv etld] (1201 Ol hereafter Paper I), where 
seven systems currently monitored by COSMOGRAIL were 
analyzed. We present here the deconvolution of a sample of an 
additional eleven lensed quasars with published time delays: 
JVAS B02 18+357, SBS 0909+532, RX J09 11.4+0551, 
FBQS J0951+2635, HE 1104-1805, PG 1115+080, 
JVAS B1422+231, SBS 1520+530, CLASS B1600+434, 
CLASS B1608+656, and HE 2149-2745. 

High-accuracy relative astrometry of lensed quasars allows 
one to test the ability of a smooth mass model of the galaxy 
to reproduce the observed image configuration. The failure of a 
smooth lens model to do this, which means that we identify an 
astrometric anomaly, may be caused by a perturbation of the po- 
tential by dark matter clump(s) along the line-of-sight, and likely 
located in the halo of the main lens. During the past five years, 
the occurrence and amplitude of astrometric anomalies caused 
by dark matter clumps has been questioned and their use to study 
the amount of substructure s in galaxies has been investigated 
dZackrisson & R iehm 2009, Paper I and references therein). We 
use the seven quads studied here and in Paper I plus seven other 
quads with astrometry derived with ISMCS or VLA+HST imag- 
ing, to revisit the evidence for astrometric anomalies. We also 
quantitatively address the role of the nearby lens environment in 
producing these anomalies. 

The deconvolution method we used also provides accurate 
morphological information on the lensing galaxies. With these 
data, we compared the mass and the light distributi on in these 
system s. This problem has first been approached bv lKeeton et al.l 
d 19981) who compared the mass and the light distribution for 14 
multiply -imaged quasars (seve n do ubles and seven qu ads). More 
recently, iFerreras et al.l (120081) and iTreu et all (120091) studied a 
samp le of nine (resp. 25 ) lensing galaxies from the SLACS sam- 
ple (iBolton et al.1 120061) . Although these works probe different 
populations of lenses, the "SLACS lenses" lying in typical ellip- 
tical galaxy en vironments while lensed quasars often lie in richer 
enviro nments dHuterer et al.ll2005tlOgu ri et al. 2 0051: ISuvu et al.l 
2009), they find that the PA of the light and of the mass do agree 
within 10 degrees. Conversely, no correlation has been found be- 
tween the ell i pticity of the light and of the mass distributions in 
iKeeton et al.1 dl998h and in lFerreras etai1(l2008l) . 

The outline of the paper is the following. Section |2]presents 
the eleven lensed systems deconvolved with the ISMCS tech- 
nique. The deconvolution method and the derived astrometry are 
described in Sect. [3] Section|4]explains the mass modeling strat- 
egy. We discuss the lens models for doubly imaged quasars in 
Sect. In Sect. [6] we present an enlarged sample of quadruply 
imaged systems for which we provide additional mass models 
and a comparison of the shape of the mass and of the light dis- 
tribution. Finally, our conclusions are presented in in Sect. [7] 



2. Data 

We used high-resolution images of eleven gravitationally lensed 
quasars obtained with NIC2, the Camera 2 of NICMOS, Near- 
Infrared Camera and Multi-Object Spectrometer, on board the 
Hubble Space Telescope (HST). We preferred to use NICMOS 
data rather than data obtained with the Advanced Camera 
for Surveys (ACS) for two reasons. First, most of the lensed 
quasars monitored by COSMOGRAIL have been observed with 
NICMOS but not always with the ACS. Second, the field dis- 
tortion and point spread function (PSF) variations are stronger 
with the ACS than with NICMOS, with damaging conse- 
quences for our deconvolution technique, which fails to sig- 
nificantly improve the Tiny Tim PSF for ACS images. The 
NIC2 data were obtained in the framework of the CASTLES 
project (Cfa-Arizona Space Telescope LEns SurvejQ, PL: E. 
Falco). Each data set is composed of four (eight in the case of 
JVAS B 1422+23 1) dithered frames acquired through the F160W 
filter (H-band) in the MULTIACCUM mode, each of them is 
a combination of 18 to 20 samples. Table Q] summarizes some 
useful information about the systems: the coordinates in J2000, 
the type of configuration, the redshift of the lens zi, the redshift 
of the source z s , and the angular separation (measured on the 
NICMOS images, see next section). Table [2] displays informa- 
tion about the image acquisition: the date of observation of ev- 
ery object, the number of frames NF and the total exposure time 

*exp ■ 

Additional details about the lensed quasars under investiga- 
tion are provided below. 

- JVAS B0218+357 (a). This system is a doubly imaged 
bla zar (BL Lac type), first considered as such in 1 19931 
by iPatnaik et alj in the context of the Jodrell Bank/VLA 
Astrometric Survey (JVAS). It has the smallest image 
separation found so far in l ensed quasar system s. The lens is 
a gas-rich late-type galaxy dKeeton et alj|1998l) . 

- SBS 0909+532 (b). Discovered as a quasar by 
IStepanvan et al.l d!991l) in the framework of the Second 
Byurakan Survey (SBS), the hypothesis of its le nsed nature 



was first pro posed by Kochane ket al. _l 199% and then 



supported by IOscozetal.1 d!997l) . In 119981 lEngels et all 



1 MCS is an abbreviation build based on the name of the three de- 
signers of the method, Magain, Courbin, and Sohy. 



resolved the system into a pair of point sources during the 
Hamburg-Cfa Bright Quasar Survey. The two lensed images 
are bright and close t o each other, which made it difficult 
for iLehar et al.l d2000l) to detect the lensing galaxy even on 
high-resolution near-IR images: they concluded that the lens 
is probably an early-type galaxy with a large effective radius 
and low surface brightness. 

- RX J091 1.4 +0551 (c). T his lensed quasar was discovered 
in 1 19971 by iBade et all amongst the AGN candidates of 
the ROSAT All-Sky Survey (RA SS). Its lensed nature was 
confirmed by Burud et al] d!998l) . The lens is a nearly cir- 
cular early-type galaxy, which belongs to a massive cluster, 
centered 38" south-west of the lensed system. It is also 
known to have a small sat ellite galaxy in the north-west 
direction dKneib et alJl2000h . 

- FBQS J095 1+2635 (d). This system was dis covered during 
the FI RST Bright QSO Survey (FBQS) and lSchechter et all 
d 19981) first identified the two observed objects as lensed 
images of the same source. The lens is an early-type galaxy 



http : //www . cf a . harvard . edu/castles/ 
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(f) PG 1115+080 


ll h 18 m 17 s .00 


+07°45'57'.'70 
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0.3098 ± 0.0002* 


1.722' 


2.43 


(g)JVASB 1422+231 


14 h 24 m 38 s .09 


+22°56'0O760 


Q 


0.3366 ± 0.0004"' 


3.62" 


1.28 


(h) SBS 1520+530 


15 h 21 m 44 s .83 


+52°54'48'.'60 


D 


0.761° 


1.855 ±0.002" 


1.57 


(i) CLASS B 1600+434 


16 h 01 m 40 s .45 


+43° 16' 47/80 
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0.41° 


1.59' 


1.40 


(j) CLASS B 1608+656 


16 h 09 m 13 s .96 


+65° 32' 29".' 00 
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(k) HE 2149-2745 


21 h 52 m 07 s .44 


-27 o 31'50720 
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0.603 + 0.001" 
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1.70 



Notes: References for the redshifts of the lenses and sources: (a, lBrowne et al . 1993; S tick et & Kuehd 1993: Carilli et al. 1993 (b) Cohen et alj2003l (cl lLubin et all 
120001 (d) lKochanek et alJI 19971 ; IQscoz et al]| 19971 (e) Kneib et al. 2000 (f) Bade et al. 1997 (g) Eigenbrod et al. 2007 (h) Jakobsson et al. 2005 (i) Lidma n et al] 
2000 (j) Smetteetal. 1995 (k) Tonrv 1998 (1) Wevmann et al. 1980 (m) Tonrv 1998 (n) Patnaik et al. 1992 (o) Auaer et al. 2008 (p) Chavushvan et al. 1997 
(q) Hansen & Hiorth 1997; Fassnacht & Cohen 1998 (r) Fassnacht & Cohen 1998 (s) Mye rs et alj 119951 ft, [Fassnacht etafl 119961 (u) lEigenbrod et al] 120071 
(v) Wisotzki eTall 19961 

Table 1. General information about the objects of the sample. Type D is for doubly imaged and Q for quadruply imaged system. 



(e.g. lEigenbrod et al.ll2007h . 

HE 1104- 1805 (ej. First consid ered as a good lens can- 
didate by IWisotzki et al.l (1 19931) in the framework of the 
Hamburg/ESO Survey (HES) for bright quas ars, the nature 
of this dou bly lensed quasar w as confirmed by Courbin et al. 
(119981) and lRemv et all (1 19981). who detected the lens, which 
is an early-type galaxy (ICourbin et al. 2000). 



PG 1115+080 (f). IWevmann etaP (Il980h first identified 
PG 1115+080 as a gravitationally lensed BAL (Broad 
Absorption Line) quasar in the framework of the Palomar- 
Green survey (PG). At that time they only det ected three 
point sources, while this system is indeed a quad dHege et al] 
1981) with two very close-by images. An IR Einstein ring, 
distorted image o f the quasar host ga laxy, was detect ed for 
the fir st time by llmpev et"al] d!998l) . According to iTonrvl 
(119981) . the lens, which is elliptical dYoo et all 120051) . is 
part of a group located at a redshift of 0.311 + 0.001. 
The evidence for a dark matter substruct ure in th i s lens 
was studied bv lChiba et all d2005l) . while lYoo et al] d2006l) 
investigated its halo morphology. They found an upper limit 
of the order of 1 mas on the possible offset between the 
center of the light distribution and of the halo. 

J VAS B 1 4 22+23 1 (g). Th is quadruple system was discovered 
in ll992l bv lPatnaik et alJ during the Jodrell/VLA Astrometric 
Survey. Th e lensed nature of JVA S B 1422+231 was first 
noticed by lLawrence et all d 19921) and th e n sup ported by 
lRemv etal.1 dl993l). According to ITonrvl dl998l) . as first 
publish ed by iKundic et al.l dl997l) . the lens, an early-type 
galaxy dlmpev et aljl 19961) . is part of a group at a redshift of 
z g = 0.339 + 0.002. 

SBS 1 520+530 (h). Th e Seco nd Byurakan Survey al- 
lowed IChavushvan et al.l d 19971) to discover this doubly 
imaged BAL quasar. T he lensing galaxy was identi fied by 
Cram pton et al.l d!998l) and is elliptical (e.g. lAuger et al] 
2008)'. 



CLASS Bl 600+434 (i). This syste m is a two-im a ge len sed 
quasar first identified as such by iJackson et al] d 19951) in 
the context of the Cosmic Lens All-Sky Survey (CLASS). 



Object 


Obs. date 


NF 


£ exp 


NI 




X>-,i 


(a) JVAS B02 18+357 


1997-08-19 


4 


43' 


6 


23.11 


1.00 


(b) SBS 0909+532 


1997-11-07 


4 


47' 


3 


12.70 


1.43 


(c)RX J091 1.4+0551 


1998-10-18 


4 


43' 


3 


9.47 


1.47 


(d) FBQS J095 1+2635 


1998-03-19 


4 


43' 


6 


60.37 


0.94 


(e) HE 1104-1805 


1997-11-22 


4 


43' 


4 


30.26 


1.61 


(f)PG 1115+080 


1997-11-17 


4 


43' 


5 


15.02 


2.35 


(g) JVAS B 1422+231 


1998-02-27 


8 


85' 


4 


118.3 


1.76 


(h) SBS 1520+530 


1998-07-20 


4 


47' 


5 


11.78 


1.35 


(i) CLASS B 1600+434 


1997-10-10 


4 


43' 


3 


5.54 


1.21 


(j) CLASS B 1608+656 


1997-09-29 


4 


47' 


4 


28.50 


2.09 


(k) HE 2149-2745 


1998-09-04 


4 


43' 


5 


18.49 


1.06 



Table 2. General information about the HST/NIC2 images and 
ISMCS. 



The primary lensing object is an edge-on spira l galaxy 
dJaunsen & Hiorthlll997l:lFassnacht & Cohenll 19981) . 

CLASS B1608+656 (j). This lense d quasar, which is th e 
core of a post-starb urst radio galaxy (Fassnacht et al][T996l) . 
was discovered by Mye rs et al.l d!995 ) as part of the Cosmic 
Lens All-Sky Survey. This object is quadruply lensed by two 
early-type galaxies located at the same redshift. The main 
lens belongs to a low-mass group of redshift z g = 0.631 
that is composed of eight confirmed members. Moreover, 
three other groups are located in the foreground of the 
primary lens: they contain 10 confirmed members each and 
are respect ively located at z gl - 0.265, z g2 = 0.426 and 
z g , = 0.52 (Fa ssnacht ei~alll2006l) . 



- HE 2149-274 5 (k) T his doubly lens ed BAL quasar was 
discovered in 119961 bv IWisotzki et al .1 in the context of the 
Hamburg/ESO Survey for bri ght quasars. The main deflec- 
tor remained undetec ted until | Lopez et ak| (fl 998) and is an 
ellipti cal galaxy ( e.g. lEigenbrod et aLl l2007). According to 
iFaure et al] (120041) a galaxy group lies along the line of sight 
to this system. 

Three systems with measured time delays but lying in com- 
plex or rich environm ent have been exclu ded from our analysis: 
PKS 1830-211 (e. g. lLoyell et al l Il998h. SDSS J 1004+41 12 
120071 120081) . QSO 0957+561 (e.g. 



( Fohlmei ster et alj 
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lOscoz et ai ] 1200 It [Shalyapin et al. 2008; iFadelv et all l2010h . 
Two other systems are not presented here because we already 
publis hed improved astrom etry and time delays : WFI J2033- 
4723 dVuissozet al.ll2008h and HE 0435-1223 dCourbin et alJ 

l2oToh . 



3. Iterative deconvolution 

3.1. Method 



The MCS deconvolution algorithm (Magain e tail I1998I) has 
proved to be very useful for disentangling the contribution of 
point sour ces and diffuse elements such as galaxies and arcs (see 
e.g. iBurud et al.ll2()00t iLetawe et al.l l2008). The basic principle 
of this technique is the following: it does not try to reach an 
infinite resolution in the resulting deconvolved frame, in order 
to conform to the sampling theorem. Instead, it aims at obtain- 
ing an improved resolution in the final frame in deconvolving 
not by the total point spread function (PSF) but by a partial 
kernel. Moreover, each frame is decomposed into two contri- 
butions: one from the point sources and the other one from the 
smooth structures in the form of a numerical component. In our 
case, the resolution of the final deconvolved frame, i.e. the shape 
of a point source, is chosen as a Gaussian with a full-width-at- 
half-maximum (FWHM) of two pixels, the final sampling step 
of which is two times smaller than the original one. The original 
pixel size with NICMOS depends on the observation date and 
ranges from 07075 to 07076. 

The most delicate point in deconvolution is probably the ac- 
quisition of a good PSF To ach ieve this, we used a technique 
described in iMagain et al.l d2007l) . which allows us to add a nu- 
merical component to a first approximation of the PSF so that it 
fits the actual PSF. But to do this, some PSF stars are needed. 
In our case, no star is available in the field of NIC2 which, is 
not large: 1972x1 972. That is why we used the MCS algorithm 
in a slightly diff erent way, i.e. coup l ed wi th an iterative strat- 
egy described in IChantrv & Magainl d2007l) and called ISMCS. 
Here is how it works. All available frames for one object are si- 
multaneously deconvol ved a first time with PSFs synthesized by 
the Tiny Tim software dKrist & Ho ok 2004) aimed at producing 
PSFs for every configuration of the HST These PSFs are sat- 
isfactory as a starting point, but not good enough for the level 
of accuracy we aim to reach (see IChantrv & Magain 2007} for 
more details). This first deconvolution produces a first estimate 
of the diffuse structures which, once reconvolved to the initial 
resolution, can be subtracted from the original frames. We thus 
obtain new images with point sources less contaminated by ex- 
tended components. We are then able to adjust our PSFs on these 
point sources and obtain a new set of PSFs that we once again 
use to simultaneously deconvolve our frames: we obtain a new 
background map (more accurate than the previous one) that we 
can subtract once more from the original frames. This process 
has to be repeated until a reasonable reduced chi square, is 
reached on the image area of interest. In practice, we stopped the 
process when < 0.2 between two consecutive iterations or 
when^ < 1.0. This^ is defined as follows: 



N 



' M{x)-D{x) 
o-(x) 



(1) 



where N is the number of pixels in the image or in the area of 
interest, M(x) is the solution reconvolved by the partial PSF, 
D(x) is the observed signal, and cr(x) is the standard deviation 
associated to that signal. 



3.2. Results 

The number of iterations NI needed to reach convergence with 
ISMCS and the^ after the first and last iterations (resp.^^ and 

X 2 rl ) are presented in the three last columns of Table|2] The origi- 
nal frame, i.e. the combination of the original F160W images, as 
well as the deconvolved frame, are displayed in Fig. Q] for each 
studied system. The labels of the lensed images are those of pre- 
vious studies. These results were obtained with a combination of 
a smoothed numerical background and an analytic galactic pro- 
file, which allowed us to obtain complementary constraints on 
the shape of the light distribution in the lensing galaxies, i.e. PA 
of the semi-major axis in degrees positive east of north, the el- 
lipticity e (e — 1 — bid), the effective semi-major axis (a e ffX the 
effective semi-minor axis (b e g) and the effective radius (R e ff)- 
The latter is the geometric m ean between the effective major 
and minor axes of the ellipse (Kochanek 2002). These param- 
eters are listed in Table [3] with their +lcr internal error bars, i.e. 
inherent to the deconvolution process (see below how they are 
obtained). As all the results, they come from the last iteration of 
ISMCS, but this time performed without the numerical channel 
to ensure that the maximum amount of light of the galaxy is in- 
cluded in the analytical p rofile. For early-type galaxies, we fitted 
a de Vaucouleurs model dde Vaucouleursl 19481) and for late-type 
galaxies, i.e. for JVAS B0218+357 and B 1 600+434 we fi tted an 
ellipsoidal exponential luminosity profile dFreemanl fl970). 

Also displayed in Fig. [T]are the mean residual maps^ from 
the first iteration, i.e. with Tiny Tim PSFs, and last iteration of 
ISMCS, i.e. with improved PSFs. Their color scale ranges from 
-5<x in black to +5cr in white. The black rectangle drawn on each 
map represents the zone used to estimate the;^- The orientation 
of these rectangles is different for each object as it is the one of 
the original images. The application of ISMCS removes nearly 
all artifacts left by the Tiny Tim PSFs. Though there are still 
some structures in the residual maps from the last iteration, there 
is no systematic structure underneath the point sources: another 
iteration of ISMCS would therefore not be able to improve the 
residuals, which are likely caused by PSF variations with the 
position on the detector. 

We list the astrometry and photometry of the point sources 
and lensing galaxy in the Vega system derived using our de- 
convolution technique in Table [4] These quantities are derived 
from the deconvolution of the individual frames using the PSF 
retrieved with the iterative method. This technique allows us to 
estimate the +lcr error bars on the relative astrometry and pho- 
tometry, to correct the relative positions from the known distor- 
tions of the NIC2 camera frame by frame, as well as from the 
difference of pixel scale between the x and y directions. The er- 
ror bars derived this way are internal error bars inherent to the 
images and their processing. To estimate the total error that af- 
fects our results we assumed that the main source of systematic 
error comes from residual geometric distortion. We have esti- 
mated the amplitude of the residual error caused by distortion 
(i.e. our maximum total errors MTE) by comparing the rela- 
tive astrometry obtained with various HST data of the Cloverleaf 
dChantrv & Magainll2007h . Here, we derived the MTE of a sys- 
tem by scaling the MTE found for the Cloverleaf with the image 
separation of the considered lensed system. This is a maximum 
error because it is estimated with the maximum spatial exten- 
sion of the system on the detector, regardless of the direction of 



3 The mean residual map is the mean of the differences between the 
model, i.e. the deconvolved frame reconvolved to the initial resolution, 
and the individual data frames, in units of the standard deviation. 
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(b) SBS 0909+532 



Fig. 1. HST/NIC2 original and deconvolved frames in logarithmic flux scale (top left and top right respectively), mean residual 
maps (linear flux scale, +5<x with white corresponding to > 5cr) from the first and from the last iterations of ISMCS (bottom left 
and bottom right respectively). 
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(d) FBQS J095 1+2635 



Fig. 1. continued. 
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(f) PG 1115+080 

Fig. 1. continued. 
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(g) JVASB 1422+231 




(h) SBS 1520+530 



Fig. 1. continued. 





(j) CLASS B 1608+656 



Fig. 1. continued. 
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(k) HE 2149-274 

Fig. 1. continued. 



Object Label 8 e (°) e a eff (") £> cff (") R cfi (") 



(a) JVAS B02 18+357 


G 


-73.60 


±2.20 


0.04: 


± 0.007 


0.32 


± 





.01 





.30 


± 





.01 


0.31 


± 





.01 


(b) SBS 0909+532 


G 


-48.10: 


t 16.90 


0.11 


±0.08 


0.57 


± 





.06 


0. 


,51 


± 





.04 


0.54 


± 





.02 


(c)RXJ09 11.4+0551 


G 


-70.00 


±4.80 


0.11 


±0.01 


1.08 


± 





.01 


0. 


.97 


± 





.01 


1.02 


± 





.01 


(d) FBQS J095 1+2635 


G 


12.80: 


± 2.20 


0.47 


±0.03 


1.00 


± 





.04 


0. 


60 


± 





.01 


0.78 


± 





.01 


(e) HE 1104-1805 


G 


54.70 


±7.60 


0.14 


±0.03 


1.05 


± 





.02 


0. 


.91 


± 





.02 


0.98 


± 





.01 


(f)PG 1115+080 


G 


-67.50 


±0.50 


0.04 


±0.01 


0.94 


± 





.01 


0. 


.90 


± 





.01 


0.92 


± 





.01 


(g) JVAS B 1422+231 


G 


-58.90 


±0.80 


0.39 


±0.02 


0.51 


± 





.02 


0. 


.33 


± 





.02 


0.41 


± 





.02 


(h) SBS 1520+530 


G 


-26.50 


± 1.00 


0.49 


±0.03 


0.99 


± 





.05 


0. 


.58 


± 





.01 


0.76 


± 





.02 


(i) CLASS B 1600+434 


G 


36.90 


±2.30 


0.75 


±0.08 


0.46 


± 





.01 


0. 


.18 


± 





.02 


0.29 


± 





.01 


(j) CLASS B 1608+656 


G, 


73.50: 


±0.40 


0.45 


±0.01 


1.93 


± 





.01 


1. 


,18 


± 





.01 


1.51 


± 





.01 




G 2 


-81.10 


±0.20 


0.55 


±0.01 


1.41 


± 





.01 


0. 


.76 


± 





.01 


1.04 


± 





.01 


(k) HE 2149-2745 


G 


9.80 i 


= 5.10 


0.21 


±0.02 


1.22 


± 





.03 


0. 


,99 


± 





.01 


1.10 


± 





.02 



Table 3. Measured shape parameters for the lensing galaxies. De Vaucouleurs light profiles are used except for JVAS B0218+357 
and B1600+434 for which we used an exponential light profile. 



this extension. It varies from 0.27 to 2.61 milliarcseconds (mas, 
see fifth column of Table 0]). The internal error should be used 
over the MTE if the first is larger than the latter, which is often 
the case for the lensing galaxy. Finally, note that the magnitude 
of each galaxy is measured in an aperture equal to the effective 
radius, /? e ff, given in Table [3] 

A few remarks can be made about the results: 

- JVAS B021 8+357 (a). iLehar et all d2000h derived a posi- 
tion for the lens by averaging the results they obtained with 
the same frames as ours and with archival HST/NIC1 im- 
ages. Our position agrees with theirs because of the large 
error bars associated to their result. These authors also warn 
about additional systematic errors affecting their measure- 
ments. We may similarly underestimate the error on the lens 



galaxy position. A first indication for this is that our posi- 
tions are affected by larger internal errors than for most of the 
other objects. In addition, our position is significantly differ- 
ent fr om the one deduced from the lens model of the radio 
ring (IWucknitzetalJl2004 ARA = ff.'06, ADec = Cf.'OT) 
and from the ACS observations of I York et ail d2005l) . The 
combination of the poor NICMOS sampling with the small 
separation between the galaxy and image B, i.e. less than 1 
NIC2 pixel ~ 75 mas, precludes us from deriving accurate 
positions for B and G. The ellipticity we find for the bulge 
of this late-type gala xy agrees w i th the one extracted from 
HST/ACS images by lYork et al.1 d2005l) . However, the spi- 
ral arms are not visible on NIC2 images. Additionally, our 
image separation of 333 mas is larger than the one found by 
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Object 


Label 


ARA (") 


ADec (") 


MTE (mas) 


Magnitude 


Flux ratio 


fn\ TV A <v RH9 1 So-^^ 7 
\d) J V/\o dUZIotjj / 


A 
/\ 


U. 


U. 


n 77 

u.z / 


17.46 ± 


0.01 


1 

i . 




Id 


0.3U/J ± V.WV+ 


n 1 111 _i_ n nnn/i 
U.lz/3 ± U.U0U4 


n n 
O.z/ 


16.88 ± 


0.02 


i 7n7 _i_ n nin 
l./U/ ± O.OjO 






0. tVoU ± U.UUol 


n nn/i 1 _i_ n nnn/i 
0.094 1± o.ouy4 


n n 
O.z/ 


17.97 ± 


0.06 


/ 


/u'* odo nnnn i cqo 
(b) dbd 0909+532 


A 


U. 




n nn 
0.90 


14.53 ± 


0.01 


i . 




B 


U.yooa ± 0.0004 


n /imi i n nnno 
-0.49 15 ± O.UOUo 


n nn 
0.90 


14.67 ± 


0.01 


0.503 ± U.U14 








-U.UJ J\J ± U.UUj / 


n on 
u.vu 


1 Q AA -4 


: U. 1 


/ 


(c~) RX T091 1 4+0^1 


A, 


0. 


0. 


2.51 


17.52 ± 


0.01 


I 




A , 


U.ZU1 1 3: U.UUU7 


406Q -i- 001 


7 S1 

Z.J 1 


17.58 ± 


0.01 


\j.yj\j zc u.uuo 




A 


U.Ul JO IC U.UUU7 


n QS7S + 000R 

U.7J / J It U.UUUO 


7 S1 

Z.J 1 


1 O.Z / ± 


U.Ul 


U.J1J ZLZ U.UUJ 




B 


-2 9681 + 0013 


0.7924 ± 0.0006 


2.51 


18.60 ± 


0.02 


0.372 ± 0.004 




p. 


-0 70 1Q + 00 1Q 


n S070 + OO^Q 

U.»JUZ,U i U.UUJ7 


7 S1 

Z,. J 1 


18.25 ± 


0.06 


/ 






-1.40U1 ± 0.0 JoU 


1 1 ^to _i_ n n/nn 
l.lO/o ± U.U4/0 


Z.J 1 


/ 




/ 


fH 1 ! FRO*; TOQS1a.7^S 
l^Q^ rDyj JUVJItZOJJ 


A 
r\ 


u. 


n 

V/. 


n qo 
u.yu 


15.56 ± 


0.01 


1 

1 . 




b 


n _i_ n nm i 
O.o9o5 ± U.UUl 1 


n ciic* _i_ n nm i 
-0.O3JO ± U.U01J 


n nn 
0.90 


16.93 ± 


0.02 


O.zoj ± O.OOj 






U. /JZI ± U.UUl / 


-U.40Uj ± U.UOjV 


n on 
0.90 


17.36 ± 


0.01 


/ 


1 1 et i i A/i i one 
(e) 1104-1805 


A 

A 


0. 


0. 


O £. 1 

z.ol 


15.58 ± 


0.01 


i . 




b 


i nnon i n nnno 
Z.9UZ0 ± 0.0008 


1 i n nnoo 
-1.330Z ± O.OOzo 


O A 1 

z.ol 


17.05 ± 


0.02 


o.ziy ± U.U05 






Q73 1 + 009? 


n S 1 70 + 001 ? 

-U.J1ZU i U.UU1Z 


7 61 

Z.U 1 


17.79 ± 


0.03 


/ 


/ r\ r-ir ' 1 1 1 c , AO A 

(I) FLr 1 1 15+080 


c 


0. 


0. 


1.98 


17.13 ± 


0.01 


a ^rr . r\ Ann 

0.255 ± 0.008 




A 

A, 


i not i n nnno 
1.3286 ± O.OUOz 


t nQQO i n nm c 
-Z.033o ± O.OOIj 


i no 
1.98 


15.64 ± 


0.02 


t 

1. 






1 /1T70 i n AHTd 

1.4 //z ± 0.0024 


1 C7C7 i n nnnn 

-1.5757 ± 0.0009 


1.98 


16.10 ± 


0.01 


n/^CiC i nnii 
0.o5o ± 0.013 




b 


n Q/inc i n nnnc 
-0.3405 ± 0.0U05 


1 n^no i n nn 1 n 
-1.9595 ± 0.0010 


i no 
1.98 


17.60 ± 


0.01 


(\ 1 £LA i n f\f\1 

0.1o4 ± 0.003 






^8 1 ^ 4-0 OOA1 
U.jOIj ± U.UU+1 


1 'XAAl -a- 00J.0 

- 1 .JH-H-Z ± U.UUH-U 


1 QS 

1 .70 


16.97 ± 


0.01 




(g) J VAvS 13 1422+25 1 


A 


n io/-.f\ i n aha /i 

0.3860 ± 0.0004 


0.3 169 ± 0.0003 


1.05 


14.35 ± 


0.01 


1. 




b 


U. 


0. 


1.0j 


14.23 ± 


0.01 


1 1 1 £ i n nm 
1.1 lo ± 0.003 




r 1 


^AO -i- 000^ 
-U.jjOU ± U.UUUJ 


n?*!!^ j. n onos 

-U. / JIO ± U.UUUJ 


1 OS 

1 .UJ 


14.92 ± 


0.01 


n son -i- o 007 
u.jyu ± u.uuz 




L) 


U.V4/U ± U.UUUO 


-U.oUIZ ± U.UUUj 


i 

I.Oj 


18.11 ± 


0.02 


U.Ojz ± U.001 




yj 


n tqo 1 _i_ n nm~7 
U. i 52Y ± O.OUj / 


-0.0390 ± 0.00j4 


i ne 
I.Oj 


18.16 ± 


0.05 


/ 


<iR<i 1S70-I-S10 

^n^ jdj ijzu+jju 


A 
/\ 


n 
u. 


n 


1 7S 

1 .Zo 


17.21 ± 


0.01 


1 

i . 




b 


1.4Z/4 ± U.UUUo 


-U.OjZj ± U.UUU4 


l.Zo 


18.03 ± 


0.01 


0.4 / 1 ± u.uuz 






i i inc _i_ n nn 1 *^ 


n loi/i _i_ n nn 1 1 
-0. 3o34 ± 0.001 / 


l.ZO 


18.22 ± 


0.02 


/ 


(l) L-LAoo b 1600+434 


A 

A 


U. 


0. 


1 1 A 

1.14 


19.98 ± 


0.03 


t 

1 . 




B 


n Tinn . n nni/i 

0.7300 ± 0.0024 


1 1 On 1 i A AAA/1 

- 1 . 189 1 zt 0.0006 


1.14 


20.22 ± 


0.07 


a om i n nc£ 

0.897 ± 0.056 




Lr 


0.0044 ± 0.0040 


n o a a a i n nn /i i 
-0.8444 ± 0.0041 


1 1 A 

1.14 


19.48 ± 


0.14 


/ 


(J^ LLAoo DlOUo+OJO 


A 


n 
U. 


n 
U. 


1 71 
1 . / 1 


18.37 ± 


0.07 


1 .UoO ± U.UJO 




b 


0. /4o4 ± U.UOzO 


1 ncio i n nno a 
-1.9j /8 ± O.OOzo 


1. /I 


18.66 ± 


0.09 


t 

1. 




L, 


n h aqi. i n nmo 
0. /4o3 ± U.UOio 


c\ a a cz. i n nmQ 
-0.440J ± O.OOjj 


1 1 1 
1. /I 


18.90 ± 


0.07 


n occ i n no 1 
0.8jj ± 0.0Z1 




n 

1 7 


-l.llol ± U.UUZj 


-I.ZJZ/ ± U.UUl o 


1 71 


19.50 ± 


0.04 


U.'+lo ± U.UZH- 




G, 


-0.4561 ±0.0061 


-1.0647 ± 0.0037 


1.71 


16.85 ± 


0.01 


/ 




G 2 


0.2821 ±0.0015 


-0.9359 ± 0.0023 


1.71 


17.29 ± 


0.01 


/ 


(k) HE 2149-2745 


A 


0. 


0. 


1.39 


15.62 ± 


0.01 


i. 




B 


0.8901 ± 0.0007 


1.4461 ±0.0003 


1.39 


17.19 ± 


0.01 


0.234 ± 0.001 




G 


0.7198 ± 0.0036 


1.1498 ±0.0051 


1.39 


17.69 ± 


0.01 


/ 



Notes: (t) Possibly large systematic errors (see Sect l3.2l ; (tt) Position measured based on Gaussian fitting. 

Table 4. Relative positions, maximum total error ("MTE"), magnitudes and flux ratios of the lensed images and lensing galaxy (see 
Fig.fflfor the labels). Magnitudes are measured within 1 for the galaxy and are PSF magnitudes for the point-like sources. 



lYork et all d2005l) in the optical, i.e. 317 mas, but really close 
to the 334 mas found in the radio. 
- SBS 0909+532 (b). This is the first published measurement 
of the lensing galaxy posi t ion an d morphology. Contrary 
to the claim of Lehar et al] (|2000), we do not find that the 
galaxy has a low surface brightness. Moreover, its effective 
radius is fairly small, i.e. 0'.'54. Regarding the residuals, 
we notice that at a long distance from the core, the spikes 
of the PSF of the bright point sources, which are not very 
well sampled (they have a FWHM of less than two pixels), 
appear in the background and on the residual map. This is 
because (1) these spikes are not accurately modeled by the 
Tiny Tim software which we use to create a first guess of the 
PSF, (2) the lensed images of SBS 0909+532 are brighter 
than in most other systems (3) we did not attempt to correct 
the PSF at distances which are much longer than the extent 



of the whole lensed system. 

- FBQS 0951+2635 (d). The ellipticity we find with ISMCS, 
i.e. e - 0.47 + 0.03 , is higher than the one found by 
Uakobsson etal] (120051) . i.e. e = 0.25 ± 0.04. The latter is 
likely underestimated because of a too large smoothing 
parameter or an oversimplified PSF because they used the 
MCS deconvolution algorithm but not combined with an 
iterative strategy. 

- CLASS B 1600+434 (i). The high blending between the lens 
and image B prevents us from reaching the desired level of 
accuracy on the point source position, and the error bars 
inherent to the deconvolution are larger than the estimated 
maximum total error (MTE) for the right ascension of image 
B. 
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Additionally, the deconvolution is unsatisfactory, although 
better than with a pure Tiny Tim fit, for three systems, i.e. 
HE 1104-1805, PG 1115+080, and CLASS B1608+656. The 
source in CLASS B1608+656 might not be exactly point-like 
in the near-IR, as suggested by the remnant structures in the 
residual map and because it is not an AGN, but the core of a 
post-starburst radio galaxy. Moreover, as for PG 1 1 15+080, we 
highlight a partial Einstein ring, which likely contributes to the 
loss of accuracy for these two objects: it is very bright and can af- 
fect the PSFs. For HE 1 104-1805, the angular separation is larger 
than in most cases and, as a consequence, the effect of the distor- 
tions must be stronger. Indeed, the distortions differ according to 
the position of the point sources on the detector and to the po- 
sition of the object on the detector, which varies from image to 
image as they are dithered. This probably results in a disagree- 
ment between each frame and consequently in larger error bars 
because they come from the individual deconvolution of each 
frame. We note that the convergence criterion was respected in 
the three cases and that the last iteration did not improve the re- 
sults. 

We also note that our results for the lensed images posi- 
tion are compatible with those found in the CASTLES database 
within their error bars (our error bars are smaller by typically 
a factor > 2), which is not always the case for the position of 
the lensing galaxies. Indeed, these extended structures are more 
sensitive to the shape of the PSE 

4. Lens modeling 

We modeled the ma ss distribution of the lensing galaxies as in 
IChantrv et aLl([2010l) , using the LENSMODEL software vl.99o 
dKeetonll200iX Two different mass models were systematically 
conside red for the main lensing galaxy. First, an isothermal 
profile (Kassio la & Kovnerlll993l: lKeeton 2001a) with a surface 
density profile of the form 



b 



2(x 2 +y 2 /q 2 ) 1 ' 2 ' 



(2) 



with q being the axis ratio and b a norm alization factor. 
Second, a de Vaucouleurs profile (Keeton 2001a) of the form 



k = kq exp 



+ y 2 /g' 
R, 



2 a/4 



(3) 



where k=7. 66925 is a constant, kq is the mass scaling pa- 
rameter and R e is the effective radius. For the two systems 
JVAS B0218+357and B 1600+434, we also used an exponential 
profile (lKeetonll2001al) that has a surface density of the form 

ko 17 x 2 +y 2 /q 2 \] 

where the parameters are similar to those described above. The 
effect of the environment on the gravitational potential is ac- 
counted for using a shear term characterized by an amplitude y 
and orientation 9 y , in our conventions, 9 y points toward the mass 
producing the shear. The agreement between the model and the 
parameters is quantified through a x 1 term calculated in the im- 
age plane: 



X 



S 



(PobsJ Pmod,i) 



(5) 



where pi are all the observable quantities that are used to 
constrain the model and <x,- are the associated standard errors. 
The x 2 can furtherm ore be expressed as a sum of four contribu- 
tions (Keeton 2001b): 



2 _ 2 .2 .2 .2 

X ~ X img /^lens Xmorph Xf\ u 



(6) 



where x\ mg encodes the contribution of the astrometry of the 

lensed images and x 2 ens is the contribution of the astrometry of 
the lensing galaxy, in both cases using MTE as error bars when 
they are larger than intrinsic errors. The term^^ o h refers to the 
contribution of the morphological parameters of the lens, namely 
the ellipticity e, the position angle 9 e and the effective radius 
R e /f (Table |3). This term is different from zero only when we 
use non-spherically symmetric mass model, i.e. singular isother- 
mal ellipsoid (SIE), de Vaucouleurs (DV), and exponential-disk 
(ExpD). When we modeled quads with SIE models, we used a 
large error for the mass ellipticity (cr e — 1) to allow it to deviate 
from that of the light. We used the uncertainty associated to the 
observed light (Tab. |3]l profile in all the other cases. When we 
modeled doubles, in addition to elliptical mass distributions, we 
also considered singular isothermal sphere (SIS) models. 
The last term of Eq. [6] i.e. x\ , has the form 



2 

XfLu 



z 



(fj-MjUc) 2 



(7) 



f,j 



where the index j denotes individual lensed images that have an 
observed flux fj + cr/j. The lens model gives the magnification 
Mj = |det(jU/)| of an image j, which has a magnification tensor 
Hj, and the intrinsic flux of the source is noted f src . To avoid that 
our models for doubles became underconstrained, we used the 
flux constraint through the measured NIR flux ratio of the two 
point sources with a lcr error of 10% to ac count for micro lens- 
ing, dust extinction, and time-delay (lYonehara et al.l l2008). For 
two systems (i.e. RX J091 1.4+0551, B1608+656), a second 
lensing galaxy is visible in the images and is included in the 
model as an SIS, assumed to be at the redshift of the main lens- 
ing galaxy. The position of this second deflector is fitted to the 
observed position using the uncertainties displayed in Table [4] 
All models were computed for a flat universe with the following 
cosmological parameters: Hq = 70 km/s/Mpc, Q. m = 0.3, and 
Q A = 0.7. 

The strategy used to search for the best model and to estimate 
the uncertainties is the same as in Paper I. First local and global 
minima were searched for by running the minimization from an 
ensemble of at least 1000 randomly distributed initial param- 
eters. Second, we sampled the parameter space in the vicinity 
of the identified global minimum using an adaptive Metropolis 
Hastings Monte Carlo Mar kov Chain (MCMC ) algorithm imple- 
mented in LENSMODEL dFadelv et al.ll2010l) . For each point of 
the MCMC, we calculated the relative likelihood of a parame- 
ter p based on the^- 2 statistics (i.e. L(D\p) = exp(-y 2 /2)), and 
calculated a 68% confidence interval for each parameter. 

The parameters of the best-fit models are displayed in Table 
[5] The columns display the following items: the name of the ob- 
ject, the type of mass distribution used, the mass scale parameter 
(i.e. the angular Einstein radius in arcseconds), the mass 
distribution ellipticity e and its orientation 9 e in degrees positive 
east of north, the effective radius R e fj in arcseconds for a DV or 
an ExpD model, the intensity of the shear y and its orientation 
9 Y in degrees (east of north), the model flux ratios, the number of 
degree(s) of freedom (d.o.f.), the^- 2 of the fit, and the predicted 



D. Sluse et al.: Deconvolution and modeling of 25 lensed quasars 



13 



time delays in days. For the quads and in the same column as 
the^- 2 , we also give thexf ma (resp. x^ ens ) as defined in Eq.|6] A 
positive time delay, Af^B > 0, means that the flux of A varies 
sooner than that of B. The median value of each parameter and 
its 68% confidence level is shown in Table [6] 



5. Doubly imaged systems 

Owing to the limited number of constraints provided by doubles, 
their use to derive Hq is still a matter of debate. On the one hand, 
owing to the generally large opening angle between the lens and 
the two images (i.e. opening angle close to 180° when the lens 
and the images are nearly aligned), the time delays depend little 
on the structure of the quadrupole ter m of the potenti al, and so 
on the intrinsic ellipticity of the lens (Kochanek 2002). On the 
other hand, the convergence produced by the lens environment 
may lead to large systematic errors on the time delays, which 
may be difficult to identify and disentangle from intrinsic ellip- 
ticity because the few number o f constraints provided by dou- 
bles (Kee ton & Zabludofnl2004l) . Although it is generally true 
that intrinsic and extrinsic sources of shear are difficult to dis- 
entangle in those systems, this is not always the case. From the 
comparison of an SIS+y model (i.e. spherically symmetric po- 
tential for the lensing galaxy) with a DV+y and with an SIE+y, 
which encode information about the lens morphology and al- 
low one to probe two different steepnesses of the mass distri- 
bution, we test how the intrinsic ellipticity affects the shear es- 
timate. We then check whether the position angle of the shear 
By, derived from our models, is supported by published envi- 
ronmental studies. Finally, we propose a ranking of the systems 
based on our ability to disentangle the shear caused by intrinsic 
ellipticity and nearby/line-of-sight environment. For this analy- 
sis we considered six out of the seven doubles presented above, 
excluding B0218+357 for which possible systematic errors on 
the galaxy position may strongly bias the models, and the four 
systems HE 0047-1756, SDSS Jl 155+6346, SDSS J1226-0006, 
and HS 2209+1914 presented in Paper I. A table compiling all 
our models (including the new SIE models of the four doubles 
analyzed in Paper I) for doubles is presented in AppendixlAl 

Although the DV and the SIE models encode information re- 
garding the intrinsic structure of the lens, the shear position an- 
gle predicted by these models is in general similar (within ~ 10°) 
to the PA of the shear predicted by the SIS model. Only three 
out of ten doubles of our sample (SBS 1520+530, HE 2149- 
2745, and CLASS B1600+434) show evidence for a different 
9 y once the observed ellipticity of the lens galaxy is included 
in the models. FBQ J095 1+2635 also shows a significant dis- 
agreement once the SIE model is considered. Assuming that (e, 
9 e ) provide a good proxy on the structure of the mass profile 
(see discussion in Sect. 16.31 ). the lack of change in 6 y suggests 
that the internal galaxy structure is not the main source of shear 
in most of the systems. This is also supported by the generally 
small change (typically < 15%) in the shear amplitude between 
the SIS and SIE in these systemfl The difference of shear am- 
plitude between the SIS and the D V model is well understood 
and explained in iKochanekl (11991). 

For the three systems where the internal structure has a sig- 
nificant effect on 9 r , we checked if the 9 7 provided by the ellipti- 
cal models is more realistic than the one provided by the SIS+y 
model: 



- For SBS 1520+530, 9 7 derived from the elliptical models 
points roughly toward the object M identified in Fig. Q] the 
agreement is better for th e SIE than f or the DV model. This 
object was suggested by Faure et al] d2002h to be an impor- 
tant source of shear together with a galaxy cluster/group lo- 
cated about WSW of the lens. Therefore, it sounds plausible 
that these two overdensities are the main sources of external 
shear. 

- The case of HE 2149-2745 is more subtle as the shear asso- 
ciated to the SIS model is less than 2%. This does not au- 
tomatically mean that HE 2149-2745 offers an unperturbed 
line of sight but that either several sources of shear cancel 
out or that the lens lies at the cente r of a group/cluste r, as in- 
deed suggested by other studies. iFaure et al] d2004l) . based 
on a weak-lensing analysis of the field of view, found an 
overdensity in the cluster map that they associated with a 
cluster nearly centered on the lens. The spectroscopi c study 
of the field of view carried out by Momch eva et alj (|2006) 
shows that there are in fact three major galaxy groups on the 
l.o. s. toward HE 2149-2745. These authors found that the 
total shear produced by these groups at the location of the 
lens is (y, r )=(O.O22, -68°). The good agreement with the 
value of 9 y derived with both the SIE and the DV models 
(0 r (DV) ~ -72.5 ± 5, r (SIE) ~ -82.4 ± 7) may not be 
fortuitous and indicates that for this system proper account 
of the environment and intrinsic ellipticity are primordial, a 
conclusion that was not obvious based on the results of the 
SIS+y model only. 

- CLASS B1600+434 is more difficult to assess. The amount 
and direction of the shear derived with the ExpD, DV, 
and SIE models are si gnificantly different. The environ- 
ment h as b een studied by Williams et all d2006l) . lAuger et al.l 
d2007l) and iFassnacht et aL ( 2008b based on imaging, spec- 
troscopic, and X-ray data. These studies unveiled that 
CLASS B 1600+434 is part of a small galaxy group (cr v ~ 
lOOkm/s), too faint to be detected in the X-ray and with 
only seven confirmed members (including the lens). The 
main contributor to the external shear seems to be a rela- 
tively "bright" SBa galaxy located about 4.35"SE of the lens. 
The other group members seem to produce a shear mostly 
in the direction 9 ~ 82° dAuger et al.l 120071) . which devi- 
ates from any of our estimates of 9 y . Because the spectro- 
scopic identificatio n of group candidates may be incomplete 
dAuger et al.l 120071) . the global effect of the group remains 
uncertain. Owing to the small number of group members, it 
is possible that a single unidentified member may be the key 
to solve the discrepancy. On the other hand, because we do 
not know if all group members share the same halo, because 
the ellipticity of the companion is unknown and because the 
main lens is a late-type galaxy with a strickingly high ellip- 
ticity likely to differ from the halo ellipticity, a proper esti- 
mate of the external shear may be challenging in this system. 

- FBQS 0951+2635 shows a strong decrease in the shear 
amplitude when the DV model is used but an increase 
when an SIE model is considered. This increase in the 
shear amplitude also affects the shear direction, which be- 
comes incompatible with y (SIS). The value of r (SIE) ~ 
-66.8° deviates by about 30° from the main candidate galaxy 
gr oup identified about 1 a rcmin SE from FBQS 095 1 +2635 
by IWilliams etail {2006) in the direction 9 —30°. Other 
groups identified by these authors might also play a role. 



4 HE 0047-1756, and in a less significant way SBS 0909+532, show In summary, heterogeneous environment information exists 

the strongest changes in y. for six out of the ten systems. For two of them (SBS 1520+530, 
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Notes: (t) The models of this lensed quasar might be affected by strong systematic errors owing to the large uncertainty on the exact location of the lens from the 
NICMOS data, (tt) an exponential profile better reproduces the light distribution of the spiral lens galaxy and therefore the ExpD mass model should be prefered to 
the DV mass model. The i? c ff used for our DV mass model is the effective radius derived in Tab.[3]usmg exponential light distribution. 

Table 5. Results of the parametric modeling of the 1 1 newly deconvolved quads. Quantities which are results of the modeling (i.e. 
not explicit model parameters) are displayed in italic and parameters which are fitted to their observed values are displayed into 
parentheses. Note that flux ratios are used to constrain the lens model only for doubles. 



HE 2149-2745), our simple models seem to catch the effect of 
the environment and internal structure well. For two other sys- 
tems (HE0047-1746 and J1226-006), a nearby galaxy produces 
a fair fraction of the shear needed by the model (see Paper 
I for discussion). Therefore, we think that these four systems 
are those that are better suited for time-delay studies such as 
COSMOGRAIL. For B 1600+434, despite of the low density of 
the lens environment, our simple models miss to properly iden- 
tify the source of shear. Four of the remaining systems show ei- 



ther a very high amplitude of shear without obvious contribution 
in the field (Jl 155+6346, HE 1 104-1805), a likely small contri- 
bution of the environment (HS 2209+1914) or a significant de- 
pendence of the shear with the exact mass profile of the lensing 
galaxy (FBQS 0951+2635). The last system (SBS 0909+532) 
has 6 e nearly orthogonal to 6 y such that the intrinsic and extrin- 
sic source of shear will be hard to disentangle. These last six sys- 
tems, with the exception of HS 2209+1914, which needs more 
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Notes: (t) The models of this lensed quasar might be affected by strong systematic errors owing to the large uncertainty on the exact location of the lens from the 
NICMOS data, (tt) an exponential profile better reproduces the light distribution of the spiral lens galaxy and therefore the ExpD mass model should be prefered to 
the DV mass model. The R^g used for our DV mass model is the effective radius derived in Tab.[3]using exponential light distribution. 



Table 6. Median value of the model parameters and 68% confidence interval of the 1 1 newly deconvolved quads. 



detailed observations, seem less well suited for deriving Ho from 
time-delay studies. 

6. Quadruply imaged systems 

The four quadruply imaged lensed quasars we modeled (Tables|5] 
and [6]) are poorly fitted by the DV+y model while the SIE+y 
model provides an acceptable fit to the observed image po- 
sition. Still, the observed position of the lensing galaxy in 
PG 1 1 15+080 and JVAS B1422+231 is not well reproduced by 
the SIE+y model, which requires a lens offset by resp. 0.02" 
and 0.01" (i.e. 68 and 37 hr l pc) w.r.t. the observed position. 
As pointed out earlier in the literature, it is necessary to account 
explicitly for the lens galaxy environment in these systems to re- 
produce both the lensed images and th e lens galaxy position (e.g. 
Ilmpev et al.lll998HKundic et al.ll 19971) . The modeling of the four 
quads studied above illustrates the shortcuts of single mass mod- 
els and advocates for mass models that account explicitly for the 
local environment of the lens. Therefore we decided to model 
most of the quads with milli-arcsec accuracy astrometry, taking 
into account the effect of satellites/companion galaxies detected 
in the images. We focused on the effect of the objects located 
in projection in the vicinity of the lens that contributes to the 
higher order terms of the potential, and are therefore more likely 
to influence the image positions. 



6. 1 . New models and extended sample 

We gathered a sample of 14 quads that we modeled with an SIE 
model for the main lensing galaxy, and a SIS model for any vis- 
ible companion/satellite galaxy or nearby group. The remain- 
der of the environmental perturbation to the gravitational poten- 
tial was modeled as an external shear. The modeling method- 
ology is similar to the one outlined in Sect. [4] except that we 
left (e, 9 e ) free during the fit in order to compare the mass and 
the light distribution of the lens. Our sample includes all quads 
analyzed here and those studied in Paper I (i.e. RXS J 1 131- 
1231, SDSS Jl 138+0314, WFI J2026-4536). We have added 
to this sample the seven systems B0128+437, MG04 14+05 34, 
HE0435-1223, SDSS 0924+0219, H1413+117, WFI 2033- 
4723, B2045+265. In three cases (B0128+437, J0924+0219, 
Jl 138+0314), the literature data did not allow us to identify any 
obvious secondary len^|and only SIE+y models have been used. 
The data we used and the references are provided in AppendixlBl 
We always assumed the second galaxy to lie at the redshift 
of the lens. For PG 1 1 15+080 and B 1422+23 1, the second com- 
ponent of our mass model is a galaxy group. We used t he po- 
sition of the group center derived by iMomcheva et alJ {2006) 



5 iFaure et al ] d201lh identified a bright substructure near the center of 
the lens J0924+0219, but they demonstrate this is very unlikely to be a 
nearby satellite galaxy. 
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Notes: (f) This model uses the group position of Momcheva et al. (2006) and the velocity dispersion of Wong et al. 12011); (tt) this model is identical to model #1 
because of the negligible contribution of the secondary lens; (t 1 1) this model uses the X-ray centroid of the group bvlGrant et al. <2004|) 

Table 7. Best SIE+y (Ni ens = 1) and SIE+SIS+y (Ni ens = 2) models for the sample of 14 quads (Sect. 16. U . The Einstein radius 
(^Ein) is in arcsec, the offset between the observed and model lens galaxy position AG is expressed in arcsec and in hT x pc and the 
RMS is in arcsec. The last column provides the contribution of the lens and image astrometry to the total y 1 (Eq.|6]l. Note that the 
flux ratios are not used to constrain the models. 



based on imaging data and we constrained the velocity disper- 
sion of th e SIS with t he vel ocity dispersion of the group ob- 
tained by I Wong et ail d201 lb . Owing to the large uncertainty 
on the group position, we also tested an alternative model for 
B 1422+2 31 with a group l ocated at the centroid of the X-ray 
emission (Grant et al. 2004) but keeping the velocity dispersion 
of the SIS unconstrained. For PG 1115+080, we have not im- 
plemented this alternative model because of the more diffuse 
X-ray emission, which leads to a very u ncertain group centroid 
dGrant et alJl2004tlFassnacht et alJl2008l) . 

The parameters of the best models we found are presented 
in Table [7] Models with one (Ni ens — 1) and two (Ni ens — 2) 
mass components are presented. When Ni ens — 1 (resp. Ni ens — 
2), we have 1 d.o.f. (resp. d.o.f.). The quality of the fit can 
be evaluated by means of the x 1 as defined in Eq. [6] but with 
•^morph = Aflux = ^ because of our modeling strategy. We also 
provide the root mean square (RMS) between the modeled and 
observed image position and the offset AG (in arcsec and in h~ l 
pc) between the centroid of the light and of the mass distribution. 
We now use these results to test the ability of our models to 
reproduce the quads down to milli-arcsec accuracy (Sect. I6.21 i 
and to compare the distribution of the light and of the mass in 
the main lensing galaxy (Sect. [63]). 

6.2. Astrometric anomalies 

Because the errors on the astrometry vary from system to system, 
we have decided to use quantities not weighted by astrometric er- 



ror bars to identify astrometric anomalies. As a rule of thumb, we 
considered that the model disagrees with the data when the RMS 
between the model and observed image position is RMS > 2 
mas, i.e. our largest measured uncertainty on a lensed image po- 
sition. Furthermore, we considered that the galaxy is not well fit- 
ted when AG > 40 hr x pc, which i s about tw i ce th e typical offset 
between light and mass found bv lYoo et al.l (f2006) based on the 
analysis of four lensed systems with Einstein rings. This offset 
corresponds to typically 10 mas and is on average twice larger 
than the uncertainty on the lens position (except for H1413+1 17 
and B0128+437). When a single lens deflector is considered, 
the lensed images and galaxy position are well reproduced in 
6 out of 14 systems (B0128+437, HE 0435-1223, J0924+219, 
Jl 138+0314, B1422+231 and B2045+265). For four other sys- 
tems (MG 0414+0534, PG 1 1 15+080, H1413+1 17, WFI 2026- 
4536), only the lensed image positions are well reproduced, 
while for the remaining systems, both image and galaxy posi- 
tions are discrepant. The inclusion of a second lens allows a 
nearly perfect match between the model and the observations 
except for RXS Jl 131-1231, WFI 2033-4723 and H1413+117. 
For RXS Jl 131-1231 the visible faint companion would in- 
crease the anomaly if too massive, as already pointed out in 
Claeskens et al.l (120061) . The nearest galaxy in projection is lo- 
cated 15 arcsec (75 h~ l kpc at zi en s) away from the lens and may 
be un able to explain the anomaly. Current studies of the environ- 
ment (iMomcheva et alj r2006: Wo ng et alJl201lh . although they 
identify at least a foreground galaxy group, also fail to explain 
the shear needed in the models. For H1413+117, the compan- 
ion only marginally modifies the predicted centroid of the main 
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lensing galaxy but leads to a significantly lower external shear. 
For WFI 2033-4723, the companion improves the modeling but 
not enough to lead to an acceptable fit. Other nearby compan- 
ions visible on the images might be the cause of the discrepancy 
but we have not investigated this possibility because of the small 
number of observational constraints. 

A few additional remarks should be made regarding the 
modeling of the quads. First, our best models for WFI 2026- 
4536 and B2045+265 predict a value of 6 e that is significantly 
tilted compared to the observed PA. Although we cannot ex- 
clude that the gravitational potential of the main lens is highly 
perturbed by the environment, this may also be a good indica- 
tion that additional ingredients are missing in the lens model. 
This can be an effect of dark matter substructure or of visi- 
ble perturbers not accounted for by our models. The elliptic- 
ity of the second mass component and the slope of the radial 
profile might also play a role in the model. Second, even in 
the models with two mass components, the amplitude of ex- 
ternal shear reaches y > 0.2 for 6/14 systems. In three cases 
(RX J0911.4+0551, PG 1115+080, B2045+265), the modeling 
with two lenses even leads to a stronger shear than the modeling 
with a single SIE+y. This is not alarming and reflects the fact 
that the shear is not a scalar quantity and has to be added vec- 
torially. The higher amplitude of the shear observed in our sam- 
ple of quads compared to our sample of doubles may be caused 
by the increase d probability of observing a quad in rich envi- 
ronment dHolder & Schechteril20()3t iKeeton & ZabludoflE^Ol: 
iHuterer et alJ2005l) . Nevertheless, the selection bias is too strong 
to allow us to properly investigate statistical occurrence of quads 
and doubles. Third, the formally good fit of H1413+1 17 is dom- 
inated by the large uncertainty (0.02") on the lens galaxy posi- 
tion. The model indeed predicts the lens to be offset by about 
170 hT x pc (ziens = 1) compared to the observed position. In ad- 
dition, the inclusion of the closest galaxy in the environment 
only marginally improves the model and only once the pri or 
on the shear (i.e. y < 0.25, similar to Ma cLeod et al.ll2009l) is 
used. Improving the accuracy on the astrometry of the lens by a 
factor 5 would reveal whether this system is the only cross-like 
quad showing significant astrometric anomaly maybe caused by 
the environment. Finally, we should note that two positions of 
the group have been tested for B 1422+231 but onl y the group 
positio n matching the X-ray centroid identified by iGrant et all 
(2004) leads to a successful fit. In that case, the SIS model of 
the group has a velocity dispersion <x SIS ~ 436 km/s, in excellent 
agreement with the velocity dispers ion cr v - 42 km/s of the 
group derived spectr oscopically bv IWong et al.l (1201 ll) . The X- 
ray group centroid of IGrant et al.l d2004l) howev er disagrees with 
the position obtained by Fassnacht et al] d2008l) . which is nearly 
centered on the main lens. 

In summary, astrometric anomalies are frequent when a sin- 
gle galaxy is used for the modeling. Nevertheless, in most of 
the cases, the inclusion of a nearby companion/group in the lens 
model led to a much better model. This effect of a massive com- 
panion on the image positions and the way the lensing galaxy 
could absorb t he perturbation ind uced by the latter agrees with 
the findings of IChen et al.1 (|2007) based on the study of mock 
lenses. For two systems, WFI2033-4723 and RXS Jl 131-1231, 
the companion galaxy is not sufficient to provide a perfect fit 
to the astrometry. In the case of WFI 2033-4723, other galax- 
ies projected at an angularly small distance from the lens (i.e. 
typically < 10 arcsec) might explain the observed discrepancy. 
RXS Jl 131-1231 remains the clearest case in our sample of an 
astrometric anomaly. Its amplitude is consistent with a pertur- 
bation produced by a dark matter substructure located near the 



Einstein radius of the lens (IChen et al.l l2007). Problems encoun- 
tered in reproducing the flux ratios, which are a lso affected by 
significant microlensing (Slu se et aLll2006ll2007l). and problems 
to reproduce the preliminary time delays i|Morgan et all [2006 ; 
IKeeton & Moustakasll2009l: ICongdon et aLlboiOl) strongly sup- 
port this interpretation. Astrometric anomalies possibly caused 
by substructures might also be present in other systems such 
as WFI 2026-4536 and B2045+265, for which the best sim- 
ple model predicts a strong tilt between the mass and the light 
distributions. Finally, we mention that astrometric anomalies 
may appear when higher astrometric accuracies are considered. 
Ind eed, we did not fin d any sign of anomalies for B0 128 +437 
but iBiggs et all (120041) have shown that such anomalies are ob- 
served when sub-mas accuracy on the astrometry is used. 



6.3. Position angle and ellipticity 

We now compare the ellipticity e and position angle 9 e of the 
mass distribution of the main lens with the one derived from the 
light profile. The results are displayed in Fig.|2]for the 12 quads 
for which measurements are available. This figure confirms the 
previous fin ding that the PA of the light and of the mass generally 
agree well ( IKeeton et all 1 19981: iFerreras et ail 120081: iTreu et all 
2009). A mismatch by about 10 degrees is not uncommon but 
it is unclear whether this is real, i.e. caused by triaxiality of the 
dark matter or lens environment, or caused by errors in the PA 
measurements. The systems B2045+265 and WFI 2026-4536 
show a much larger mismatch whose origin is not yet clearly 
identified. 

More surprisingly, Figure [2] shows that the mismatch be- 
tween the ellipticity of the mass and of the light distribu- 
tion is not as large as previously claimed (Kee ton et alJ il998: 
IFerreras et al.ll2008l) . For 8 out of 12 systems, the observed el- 
lipticity of the light agrees well with that of the mass. Two 
of the four outliers (RX J091 1.4+0551, SDSS Jl 138+0314, 
B1422+231, and B1608+656) show a nearby galaxy which lies 
within the Einstein radius of the lens, an effect which could 
easily bias the measurement of the light ellipticity. It is also 
likely that for these two systems and for B2045+265, the el- 
lipticity of the second lens, which lies within 1 Re of the lens 
center, introduces internal shear not accounted for by our mod- 
els or that the system is not relaxed and poorly reproduced 
by two separated potentials^. The d iffer ence between our re- 
sults and those of IKeeton et al.l dl998l) and IFerreras et all (f2008) 
may have several origins. First, these authors mixed doubles and 
quads in their analysis, while doubles do not always provide ro- 
bust information on both the shear and the ellipticity (Sect. [5}. 
Second, our ellipticity measurements have been derived nearly 
exclusively with the same method (MCS deconvolution) applied 
to NICMOS data (only MG 0414+0534 and B2045+265 have 
morphological parameters derived with another method) while 
Keeton's measurements are more heterogeneous with possibly 
more significant bias es associated to less accurate PSF models. 
Finally, the results of Fe rreras et al.l d200 8) are derived with pix- 
elated mod els. Despite the ma ny tests and successes of this tech- 
nique (e.g. ISaha et al I l2006allbt iRead et al.l 12007b IColesI [2008; 



Leierll2009l:lLeier et al. 



201 1), it is unclear whether the uniform 



prior on the ellipticity used to sample the parameter space is 
not biasing the average pixelated mass profile derived by this 
method. 



6 This might also be the case for M G 0414+0534 but th e companion 
galaxy is possibly a foreground object dCurran et alj201 il) 
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Fig. 2. Light vs mass (model) PA (left) and ellipticity (right) for 
main lens in B1608+656 (e obs =0.45±0.01, e mod =0.84+0.01) is not 
Eintein radius of the main lens are depicted with a filled triangle, 
traces the total mass. The error bars are 68% CL error bars. 

These results suggest that the dark matter and baryon distri- 
bution in the central regions (typically 5 kpc) of lens galaxies do 
have similar shapes. This may however not be a generic result 
because of possible selection biases. First, we have probed sys- 
tems with an observed ellipticity e < 0.4, such that the effect of 
a rounder dark matter halo on the total mass distribution may be 
harder to detect. Second, it may also be that there is a selection 
bias of the lensing galaxies towards systems where the observed 
corre lations are more likely to take place (cf. Mandelbaum et all 
2009). The alignement between the dark matter and the baryons 
distribution has been found in nu merical simulations for spi- 
ral galaxies (Gusta fsson et al.ll2006l) . but still has to be demon- 
strated for massive ellipticals, where significant misalignements 
have been observed aro und one virial radius dBett et al.ll2010t 
LSkibba&Maccio1l2011h . 

7. Conclusions 

We applied the iterative strategy coupled with the MCS de- 
convolution algorithm, i.e. ISMCS, to high-resolution near-IR 
images of a sample of eleven lensed quasars. The ISMCS al- 
lowed us to obtain accurate astrometry, in most cases with 1 
to 2.5 mas error bars, and shape parameters for the light dis- 
tribution of the lensing galaxy. In three cases, the deconvolu- 
tion process was not entirely satisfactory, i.e. HE 1104-1805, 
PG 1115+080 and CLASS B1608+656, but still provided bet- 
ter results than simple PSF subtraction. For one special case, 
i.e. JVAS B0218+357, the small angular separation and the poor 
resolution of NICMOS prevents us from reaching the desired ac- 
curacy. For one system, SBS 0909+532, we were able to detect 
the lensing galaxy for the first time, which is angularly small 
and relatively bright. Two objects reveal partial Einstein rings 
(PG 11 15+080, CLASS B 1608+656). 

We provided, for these systems, simple mass models 
(isothermal and de Vaucouleurs mass distributions) that give 
an overview of the ability of smooth mass models to repro- 
duce the known lensed quasars. These models give some hints 
about the effect of the radial profile on the flux ratios and on the 
time delays. For the sample of ten doubly imaged quasars pre- 




Light Ellipticity 

our sample of quadruply imaged quasars. The ellipticity of the 
shown to ease legibility. Systems having a companion within the 
The red solid line depicts the situation where the light perfectly 



sented here and in Paper I (i.e. lChantrv et al.l2010l) . we have also 
studied the role of intrinsic ellipticity and external shear in the 
models. For half of these systemfQ (HE0047-1746, J1226-006, 
SBS 1520+530, HE 2149-2745, and HS 2209+1914) one should 
be able to unfold intrinsic and extrinsic shear and therefore min- 
imize systematic biases on Ho based on the time-delay method. 
This identification of the most promising systems is partly sub- 
jective because detailed environmental studies of the whole sam- 
ple are lacking. Ex haustive studies of lens environment are for- 
tunately ongoing dWilliams et al.l 120061: iMomcheva et al.l l2006; 
Won g et aDl201 ll) and should allow one to address this question 
more quantitatively. 

In addition to the ten doubly imaged quasars, we have per- 
formed modeling with singular isothermal ellipsoid (SIE) of a 
sample of 14 quadruply imaged systems for which we have rel- 
ative astrometry with accuracy typically better than 2 mas. This 
sample allowed us to investigate two questions: 

1. How efficient are SIE models in reproducing the image 
and galaxy position of the observed quads and how does 
the local environment of the lens influence this result ? 
We found that the environment of most of the quads is 
rich and that at least the nearest galaxy or a nearby group 
has to be included explicitly in the model to recover the 
image and galaxy positions. Following this procedure leads 
to a formally perfect fit with an SIE model of the main 
lens for 11 of the 14 systems. When this perturber is not 
included, models generally tend to predict the position of 
the main lensing galaxy significantly offset with respect 
to the observed position, the lensed image positions being 
less affected. For the two systems RXS Jl 131-1231 and 
WFI 2033-4723, we were not able to obtain an acceptable 
fit. The problems encountered with WFI 2033-4723 might 
be caused by other nearby galaxies not included explicitly 
in our model. RXS J 1 1 3 1 - 123 1 remains the best candidate 
of our sample for which astrometric anomalies may be 
caused by a dark matter substructure. Two other systems 

7 HS 2209+1914 is more uncertain because of the lack of detailed 
environmental studies 
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(WFI 2026-4536 and B2045+265) show a significant mis- 
match between the 9 e of the light and of the mass. The origin 
of this discrepancy is unclear and might also be the signature 
of dark matter substructures. We mention that astrometric 
anomalies may become proeminent when higher accuracy 
images are available. This is the case for B0128+437, which 
is perfectly reproduced by our model but shows astrometric 
anomalies once sub-mas astro metry of the lensed images is 
considered (Bi ggs et al.ll2004h . 

2. How do the morphology of the mass and of the light profile 
compare ? Our analysis of 12 out of 14 quads with measured 
structural parameters allowed us to find that the PA of 
the mass and of the light distribution are aligned within 
typically 10 degrees, in agreement with previous studie s 
<lKeeton etal]ll998t iFerreras et ail l2008t iTreu et all [2009). 
We also discovered that the ellipticity of the light generally 
agrees well with that of the total mass distribution (8 out of 
12 lenses show a good correlation). Overall, this indicates 
that the light provides a good proxy to the total mass 
distribution in these lensing galaxiefl, a result particularly 
interesting in the light of the rich environment of the 
lenses. Understanding why these correlations take place is 
important for both our comprehension of dark matter and 
galaxy formation and for time-delay studies. Therefore, an 
increase in the sample of lenses with shape measurements 
and a comparison with numerical simula tions, including 
probl e ms caused by selection biases (Ivan de Ven et alJ 
2009; MandelbaumetaLl E009), are necessary. This will 
provide important clues on the role of dark matter in shaping 
galaxies. 
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Appendix A: Models of the 11 douby imaged 
quasars 

We provide in Table IA.1I a summary of the median model pa- 
rameters and 68 % CL for the 1 1 doubly imaged quasars ana- 
lyzed in this paper and in Paper I. The SIE models for the lensed 
quasars HE 0047-1756, SDSS 1 155+635, SDSS 1226-0006 and 
HS 2209+1914 were not presented in Paper I and are discussed 
in Sect. 

Appendix B: Data for the quadruply imaged 
quasars 

We provide the information associated to the 14 quadruply im- 
aged quasars studied in Sect. [6] Table IB. II displays the relative 
RA-Dec of the individual images with their associated errors. 
Table IB ,2l provides the position, ellipticity and position angle of 
the lensing galaxy and the associated references. 
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Notes: (t) The models of this lensed quasai' might be affected by strong systematic errors owing to the large uncertainty on the exact location of the lens from the 
NICMOS data, (ft) An exponential profile better reproduces the light distribution and therefore the ExpD mass model should be prefered to the DV mass model. 
The R e ff used for our DV mass model is the effective radius derived in Tab. fusing exponential light distribution. 



Table A.l. Median value of the model parameters and 68% confidence interval for the 1 1 doubly imaged quasars. 
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0.0020 


-2.1946 


1.2601 


0.0020 


-1.4809 


1.3756 


0.0020 


-2.1128 


-0.2778 


0.0020 


8 


B2045 


0,0 


5.0e-4 


-0.1316 


-0.2448 


6.0e-4 


-0.2869 


-0.7885 


5.0e-4 


1.6268 


-1.0064 


0.0013 
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Notes : 1: This paper, 2: IChantrv etaLl )201C1), 3 : iBiggs et all 120041) . 4: iRos et all 120001) , 5: ICourbin eFai] <2010|) . 6: iKeeton etall 120061) . 7: IChantrv & Magainl 
120071) . 8: IVuissoz etail 120081) . 9: lMcKean etafl 120071) : f: 07^ =0.0038, a DeCi =0.0025; ft: or4 4 =0.0018, a DeCi =0.0018. 

Table B.l. Lensed image positions (1, 2, 3, 4) of the 14 quads studied in Sect. [6] We have shortened the names of the systems to 
ease legibility. 
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Name 


RA 01 




Decgi 


O" DeCQ\ 




°~RA G2 


Dec<32 


0" Dec G1 


e 




DA 

PA 


CTPA 


Ket 


BO 128 


0.217 


0.01 


-0.104 


0.01 


999 


999 


999 


999 


999 


999 


999 


999 


10 


TV A f ^ (\ 111 

MG0414 


1.072 


0.011 


0.665 


0.011 


1.457 


0.013 


2.122 


0.013 


0.2 


0.02 


84.0 


4.0 


4, 11 


T T I~~V \ (Of 

HE0435 


1.1706 


0.0030 


-0.5665 


0.0020 


3.7558 


f\ AAOA 

0.0020 


-4.1733 


0.0020 


0.09 


0.01 


-1.4 


1.7 


5 


RXJ091 1 


-0.7019 


0.00251 


0.502 


0.0039 


-1.4601 


0.036 


1.1678 


0.047 


0.11 


0.01 


-70.0 


4.8 


1 


J0924 


-0.1804 


0.0039 


-0.8685 


r» AAn r» 

0.0039 


999 


999 


999 


999 


0.12 


0.01 


-61.3 


0.5 


12 


PGl 1 15 


0.3813 


0.0041 


-1.3442 


f\ f\f\ A f\ 

0.0040 


3.0 


1 1.0 


-21.7 


9.0 


0.04 


0.01 


-67.5 


0.5 


1, 17 


T"» ~\7 T 1 1 1 1 

RXJH31 


-2.0269 


0.00264 


0.6095 


0.00264 


-1.93 


0.0050 


1.137 


0.0050 


0.25 


0.04 


-71.4 


2.4 


2, 14 


JH38 


-0.4633 


0.0071 


0.534 


0.0036 


999 


999 


999 


999 


0.16 


0.02 


-57.3 


6.5 


2 


in it i j 


U. 1 J>UJ 




fl ^887 


n CO 


1 87 

1.0/ 


n i 

U. 1 


A 14. 

4.11 


1 

\t. i 


QQQ 

yyy 


QQQ 


QQQ 

yyy 


QQQ 

yyy 


7 1 (\ 

1 , 1 u 


B1422 


0.7321 


0.0037 


-0.639 


0.0054 


7.0 


1.5 


-9.0 


1.5 


0.39 


0.02 


-58.9 


0.8 


1, 13 


B1608 


-0.4561 


0.0061 


-1.0647 


0.0037 


0.2821 


0.00171 


-0.9359 


0.0023 


0.45 


0.01 


73.5 


0.4 


1 


WFI2026 


-0.0479 


0.0015 


-0.7916 


0.0015 


-7.398 


0.026 


1.94 


0.0060 


0.24 


0.03 


60.8 


5.4 


2, 15 


WFI2033 


-1.4388 


0.0020 


0.3113 


0.0020 


-5.41 


0.0020 


0.285 


0.0020 


0.18 


0.03 


27.8 


4.3 


8 


B2045 


1.1084 


0.0011 


-0.8065 


0.0011 


0.4498 


0.0021 


-0.6425 


0.0021 


0.06 


0.01 


-60.4 


4.8 
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Notes: 1: This paper, 2: IChantrv et aT] 120101) . 3: iBiggs et all 120041) . 4: iRos et all 120001) , 5: ICourbin et al] <2010l) . 6: iKeeton et afl 120061) . 7: IChantrv & M again 
120071). 8:IVuissoz et al] 120081), 9:|M cKean et al. (2007), 10: Lagattuta et al. 12010), 11: Falco et al. ( 1997), 12: lEigenbrod et alj 12006]) . 13: iGrant et al] 120041) 14: 
Claes kens etal] 420061) . 15:|Morgan et al. 12004), 16: MacLeod et al. (2001), 17: lMomche va et al. 1 200^) 



Table B.2. Position of the main lensing galaxy (Gl) and of the nearest companion/group (G2) relative to the reference lensed image 
identified in Table IB. fl The observed ellipticity and position angle (E of N) of Gl are also provided. Unknown quantities have an 
assigned value 999. 



